Formation of density waves via interface conversion of ballistic and diffusive motion 
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We develop a mechanism for the controlled conversion of baUistic to diffusive motion and vice 
versa. This process takes place at the interfaces of domains with different time-dependent forces in 
lattices of laterally oscillating barrier potentials. As a consequence long-time transient oscillations of 
the particle density are formed which can be converted to permanent density waves by an appropriate 
tuning of the driving forces. The proposed mechanism opens the perspective of an engineering of 
the nonequilibrium dynamics of particles in inhomogeneously driven lattices. 
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INTRODUCTION 



Time-dependent forces are often the origin for the ap- 
pearance of complex dynamics in physical systems 
Spectacular effects such as the directed transport of par- 
ticles in the absence of a net force can be evoked by 
breaking certain space-time symmetries (parity and time- 
reversal [2]) homogeneously, i.e. through the application 
of a spatially independent driving force. In the literature 
this phenomenon is usually referred to as the ratchet ef- 
fect [3| and has attracted a great deal of attention since 
it is the working principle of molecular 0, Q or quantum 
^ motors. In the case of Hamiltonian (or deterministic) 
ratchets [7.-^] both setups with fully chaotic and 
mixed phase space 0, [l3l - [2C)| have been explored. Ex- 
perimentally, quantum ratchets have been demonstrated 
e.g. in Josephson junction arrays pH and for cold atoms 
loaded into periodically amplitude-modulated optical lat- 
tices 
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In view of the above it is an intriguing perspective to 
consider the consequences of a breaking of the parity and 
time-reversal symmetry inhomogeneously by applying a 
spatially varying driving force. Recently [23| a phase- 
modulated driven lattice with unit cells consisting of a 
few sites has been shown to provide a mechanism which 
leads to a selective patterned trapping of particles on cer- 
tain sites. The latter could be useful for atom lithography 
or as an ingredient for quantum information processing. 
Here we explore a different driven lattice reminiscent of 
superlattices in the physics of semiconductor heterostruc- 
tures but now with a spatiotemporal embedding. It con- 
sists of a lattice of domains, each containing many unit 
cells, and subject to a different time-periodic force. The 
driving is chosen such that neighboring domains possess 
oppositely directed currents ("local" ratchets). Our focus 
is the analysis of the nonequilibrium dynamics in such bi- 
modally driven devices. At the interfaces of the domains 
a controlled conversion of the dynamics from diffusive to 
ballistic behavior and vice versa is possible and can be 
understood in terms of the overlap of the Poincare sur- 
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FIG. 1. Schematic illustration of the driven lattice. Interfaces 
between the domains with different driving laws are at the 
positions {0,±NL,±2NL}. Arrows indicate the direction of 
the local currents. 



faces of section (PSS) of the individual domains. As a 
result we observe the formation of transient spatial den- 
sity oscillations. The latter persist in time if the lattice is 
fed with a constant current of particles. By switching ap- 
propriately the applied forces it is possible to manipulate 
these density waves. 



DRIVEN LATTICE HAMILTONIAN 

We consider noninteracting particles in a spatiotempo- 
rally driven lattice with potential barriers of equal height 
Vq and width I described by the Hamiltonian 
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where m is the mass of the particles and xo^i = i L 
is the equilibrium position of the ith barrier. fi(t) is 
a periodic driving law, which depends on the site in- 
dex i according to f,{t) = {-1)1^/^1 g{t) with g{t) = 
C {cos (a; t) ~\- sin (2a; t)} being a biharmonic function and 
[a;J is the largest integer k with k < x (see Fig. [1]). The 
lattice is composed of domains each consisting of N later- 
ally oscillating barriers obeying the same driving law and 
which are connected at corresponding interfaces. Each 



domain possesses the length Lb = N ■ L {N = 10^) and, 
as shown in Fig. [l] the driving law fi(t) alternates from 
one domain to the next. The amplitude C of the driv- 
ing is fixed such that — 1 < g{t) < 1, i.e. C ~ 0.57 and 
without loss of generality we assume cj = 1 and m = 1, 
as well as I = 1, Vq = 2.2 and L = 5. 

TRANSPORT PROPERTIES OF THE LATTICE 

Applying the driving law g{t) to the entire lattice leads 
to directed currents since the relevant spatiotemporal 
symmetries derived in 0] are broken. By integrating 
over the main chaotic sea with respect to p = the 
transport velocity can be estimated to vt = —0.2649 
po| . In a domain with the driving law fi{t) = g{t) par- 
ticles in the chaotic sea are on average transported in 
negative a;-direction. Since the driving changes its sign 
from one domain to the next one, i.e. fi{t) = — /i+Ar(i), 
the direction of the local current alternates, too, and 
transport is reversed. There are two types of interfaces 
that repeat in an alternating manner in the lattice. One 
type is characterized by incoming 'convergent' currents 
from both neighboring domains whereas the second one 
has 'divergent' outgoing currents pointing towards the 
neighboring interfaces (see Fig. [1]). Naturally, one could 
therefore expect that interfaces with incoming currents 
show particle accumulation whereas those with outgoing 
currents exhibit particle depletion. Although each do- 
main shows directed transport and an according drift of 
the particles, the complete lattice consisting of infinitely 
many domains does not exhibit directed currents since 
the equations of motion are invariant under the transfor- 
mation T : X ^ —X + NL, t — >■ f Yet, an ensemble of 
particles reaches this zero current limit of the complete 
driven lattice only for extremely long time scales by far 
not reached in the present work. 

RESULTS 

We explore the dynamics of an ensemble of 10^ parti- 
cles which are initially uniformly distributed in the lat- 
tice occupying 20 domains centered at the origin and 
which homogeneously fill the chaotic sea at low velocities 

v\ < 0.1. (For details of the numerical method see Refs. 

20, |2J]). The spatiotemporal evolution of the particle 
density profile (see Fig. [2]) shows due to the diffusion of 
the chaotic trajectories a significant spatial broadening. 
More importantly, we observe the emergence of a peri- 
odic stripe-like modulation of the density, i.e. for a fixed 
time the spatial distribution of the particles has alternat- 
ing flat maxima and minima corresponding to the length 
of one domain. This behavior persists until t « 5 • 10^ 
and gradually disappears thereafter. Let us analyze this 
in some more detail. 
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FIG. 2. (Color online) Particle density p{x,t). The initial 
ensemble is distributed according to p{x, t — 0) — 10"" ■ 
e{10^L - \x\) and p{v,t = 0) = 5 • 0(0.1 - (V;, = 2.2). 
For t > 10* a periodic modulation of the spatial distribution 
of the particles in the lattice is clearly visible. 



For t = 10^ (Fig- El (a)) the envelope behavior of p{x) 
is still similar to the initial distribution po{x) but there 
occur sharp dips at the positions of the interfaces. The 
dips are significantly more pronounced at the interfaces 
with outgoing local currents (inset of Fig. |3] (a)). At 
t — 10^ ( Fig. [3] (b)) the envelope of p{x) has broadened 
and the density modulations are much more pronounced 
extending over an entire domain of the lattice (inset of 
Fig. [3](b)). For x > p{x) possesses plateau-like maxima 
in domains with wt > and corresponding minima in 
domains with wt < and vice versa for a; < 0. 

To understand the mechanism which is responsible for 
the formation of the density modulation, it is instructive 
to look at the trajectories of individual particles. Fig. 
m (inset) shows the typical dynamics which consists of 
alternating phases of chaotic and ballistic behavior. Ini- 
tially, the particle propagates ballistically opposite to the 
direction of the local particle current until it reaches the 
next interface where a conversion to chaotic dynamics 
occurs. In the next domain the trajectory diffuses chaot- 
ically with a net drift according to the directed transport 
in this domain. When it arrives at the second interface a 
conversion back to ballistic behavior is observed. The dy- 
namical conversion process repeats at each interface until 
the particle starts a very long flight which lasts over sev- 
eral domains. Typically, we use the term "ballistic flight" 
for time intervals during which the particles travel over 
a long distance x > lOL in the same direction. 

To understand the interplay between ballistic and dif- 
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FIG. 3. Spatial distribution of the particles p{x) for t = 10* 
(a) andt = 10'' (b). The arrows in the maenifications indicate 
the direction of the local current in the domain. 



fusive motion at the interfaces it is useful to visualize 
the dynamics by means of stroboscopic Poincare sur- 
face of sections (PSS). A super PSS containing the de- 
tailed dynamics for times up to t 3> 10^ over the en- 
tire lattice possesses an enormously rich structure at any 
scale and is in practice inaccessible. For the descrip- 
tion of the transient dynamics characteristic for times 
t < 10^ the PSS of a domain (PSSD) is more adequate 
(see Fig. |4]). For momenta \p\ > 4.3 the dynamics is 
mostly regular, i.e. in this region the phase space is 
predominantly foliated by deformed tori of the unper- 
turbed system. A big elliptic island centered around a 
period one orbit at x ~ 0,p ~ —4.9 is encountered. Fur- 
thermore, we observe a chaotic sea confined by invariant 
curves to 3.9 ^ p ^ 4.3. Trajectories which start in this 
sea traverse the lattice domain in the direction of their 
initial velocity. Particles starting in the large chaotic 
sea (—4.3 < p < 3.9) wander diffusively through their 
domain whereas trajectories belonging to the embedded 
regular islands travel only in the direction of their initial 
momentum. Domains with driving laws g{t) respectively 
—g(t) possess mirror symmetric PSSDs with respect to 
p = 0. At the interfaces between the regions with dif- 
ferent driving laws the PSSDs of the two domains are 
concatenated. By crossing such an interface it is possible 
that a trajectory is injected from the chaotic sea into a 
ballistic island. Then the particle stays in this island and 
traverses the subsequent domain until it reaches the next 
interface where a conversion back to the chaotic sea can 
occur. Trajectories (inset of Fig. |4]) thus exhibit alter- 
nating ballistic and diffusive phases each on the scale of 
a single domain. 

Let us now return to the discussion of the evolution 
of p{x). In the beginning {t = lO**) we observe a sharp 
depletion of the particle density at the positions of the 
interfaces (Fig. [3] (a)). Every time a chaotic trajectory 
crosses an interface it is possible that a ballistic flight 
starts at this position. By analyzing the spatial distribu- 
tion of the starting positions of ballistic flights p(a;start) 
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FIG. 4. Stroboscopic Poincare surface of section 
(PSS) for a single domain of the lattice with g{t) = 
C{coa{ujt)+sm{2ujt)} (C ^ 0.57, Vb = 2.2, 1 = 1 , L = 5). 
Inset: Individual trajectory. Dotted lines: Positions of the 
interfaces with small arrows indicating the local currents in 
the corresponding domain. 



we flnd narrow peaks localized at the interfaces. Appar- 
ently, this process leads to a fast local reduction of the 
spatial distribution of the particles p{x) at these posi- 
tions since compared to diffusive trajectories, the ballis- 
tic particles travel with a high velocity away from the 
interface. However, if the local currents are incoming 
this effect competes with the accumulation of particles 
obeying chaotic dynamics. Therefore, the dip is more 
pronounced for interfaces with outgoing local currents. 

Let us now develop an understanding of the longer- 
time transient behavior leading to the plateau-like den- 
sity modulations (Fig [3] (b)). The initial ensemble is lo- 
calized in the chaotic sea of the phase space of the lattice. 
In every domain the particles experience the correspond- 
ing drift velocity and accumulate at the interfaces I with 
incoming directed currents. Even if a 'deeper' reinjection 
from the interface into an adjacent domain again occurs 
(which is a rare event) it will return to this interface 
at a somewhat later time. Hence, chaotic trajectories 
are transported several times back to the interface until 
a conversion from chaotic to ballistic dynamics occurs. 
Accordingly, diffusive particles are localized in a small 
spatial region around these interfaces I until their posi- 
tion in phase space coincides with a regular structure in 
the phase space of the adjacent domain which leads to 
a conversion into regular ballistic motion. When such 
a ballistic particle reaches the next interface II (outgo- 
ing directed currents) the trajectory can either remain 
ballistic and traverse the subsequent domain in a com- 
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paratively short time or the particle is injected from the 
ballistic island back into the chaotic sea. It is however 
important to keep track of the fact that the ballistic tra- 
jectories cross the interface II and typically enter up to 
approximately 10'^ unit cells into the next domain where 
they acquire a diffusive character. In the latter domain 
the particles experience a drift to the adjacent interface 

I, since the local current points away from the interface 

II. Obviously, this process does not occur at the interface 
I with exclusively incoming currents. Note that domains 
with left pointing current can be traversed to the right 
only by ballistic trajectories. As a result an enhanced 
number of diffusive trajectories occurs in domains with 
right pointing current for a; > and vice versa. Since the 
dwell time of chaotic trajectories in a domain is approxi- 
mately one order of magnitude larger than the dwell time 
of ballistically travelling particles this results in a peri- 
odic modulation of p{x). It is crucial to note that the 
Hamiltonian itself does not reflect such an asymmetry 
with respect to the domains with left and right-pointing 
currents corresponding to driving laws g{t) respectively 
—g{t). Indeed, the asymptotic steady state of the time- 
evolved ensemble (t ^ 10^) possesses a uniform density. 

The above nonequilibrium dynamics can be described 
by a Markov chain that represents an "interface map- 
ping" comprising the dynamics within a domain. Com- 
ponents of the Markov state are the dynamical character 
of a trajectory (chaotic / ballistic) , the direction of mo- 
tion (left / right), the position of the particle and the 
point in time it reaches the current interface. The condi- 
tional transition probabilities pij of a dynamical conver- 
sion at an interface, which constitute the transition ma- 
trices can be obtained by calculating the overlap of the 
different regular and chaotic structures in phase space 
and augmenting it by local characteristics of the differ- 
ent interfaces. The inset of Fig. [5] (compare to Fig. [3] 
(b)) demonstrates that the Markov chain indeed repro- 
duces the long-term transient spatial distribution of the 
particles in the lattice. 

It is possible to manipulate the above density profile 
such that it becomes a steady state of the spatiotem- 
porally driven lattice. To this end, we consider a finite 
open system composed of 10 domains with an incoming 
constant current of particles. After a transient phase the 
spatial distribution of the particles approaches a steady 
state (Fig. [5]) with an appearance similar to the case 
shown in Fig. |3] (b). It is now possible to render this 
steady state into a time propagating density wave by gen- 
erating a movement of the interfaces through the lattice. 
The latter is achieved by a consecutive switching of the 
driving laws of the barriers. We then observe that the 
spatial distribution of the particles follows this change of 
the lattice adiabatically, i.e. the appearance of p{x, t) is 
preserved during this process resulting in a propagating 
density wave. 




FIG. 5. Snapshots of the density p{x) in a finite, open system 
with a constant current of particles entering the lattice. Dot- 
ted lines: Positions of the interfaces at t = 10*' and t = 1.5-10® 
are identical. Inset: Spatial distribution of particles obtained 
by modelling via the Markov chain. 



CONCLUSION 

We have demonstrated novel nonequilibrium phenom- 
ena in ID driven lattices associated with the inhomoge- 
neous breaking of parity and time-reversal symmetries. 
Dynamical conversions of trajectories from ballistic to 
diffusive motion are observed at the interfaces between 
differently transporting domains in the lattice. For a 
beam of particles entering the device this mechanism 
leads to a density modulation which can be manipulated 
by switching the driving laws of the individual barriers. 
A suitable experimental setup which could demonstrate 
such density waves are multiple layered semiconductor 
heterostructures driven according to locally different AC- 
voltages or applied laser fields. Cold atoms in optical lat- 
tices formed by counterpropagating laser beams are an 
alternative experimental setup to verify our findings. By 
mounting the mirrors on piezoelectric crystals a laterally 
oscillating one-dimensional potential can be created. 
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